Parallel Information Transfer in a Multi-Node Quantum Information Processor 
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We describe a method for coupling disjoint quantum bits (qubits) in different local processing 
nodes of a distributed node quantum information processor. An effective channel for information 
transfer between nodes is obtained by moving the system into an interaction frame where all pairs 
of cross-node qubits are effectively coupled via an exchange interaction between actuator elements 
of each node. All control is achieved via actuator-only modulation, leading to fast implementations 
of a universal set of internode quantum gates. The method is expected to be nearly independent of 
actuator decoherence and may be made insensitive to experimental variations of system parameters 
by appropriate design of control sequences. We show, in particular, how the induced cross-node 
coupling channel may be used to swap the complete quantum states of the local processors in 
parallel. 



A distributed, multi-node structure has been suggested 
as a convenient means of arranging qubits in an experi- 
mentally realizable quantum computer architecture [l48| . 
Such a structure requires the ability to define an array of 
disjoint quantum processors that may be controlled lo- 
cally, with communication between local processors pro- 
vided by a coupling between nodes that may be turned 
on or off in turn. One method of satisfying these re- 
quirements is to append a small number of qubits to each 
node of an array of nearest-neighbor coupled actuator ele- 
ments. The actuator elements provide local control of the 
surrounding processor qubits and a means of transferring 
information between nodes [3]. In this work, we present 
a method for generating a universal set of gates between 
any pair of disjoint cross-node qubits via isotropic actu- 
ator couplings. 

An effective cross-node processor coupling network is 
created by taking advantage of four-body coupling terms 
between actuator and processor elements that appear in 
a manifold of excited states unused for quantum infor- 
mation storage. By moving into an appropriate interac- 
tion frame, the four-body coupling terms appear as two- 
body couplings between every pair of cross-node proces- 
sor qubits in a properly defined computational manifold. 
While this complete cross-node coupling network allows 
for a computationally universal set of operations between 
nodes, we present an explicit implementation of a paral- 
lel swap of the complete quantum mechanical states of 
two local quantum processors. Consideration of this rep- 
resentative entangling operation serves to motivate the 
broader applicability of the induced information trans- 
fer channel. Additionally, since information is never ex- 
plicitly stored for an appreciable amount of time on the 
actuators - which are exposed to higher levels of noise 
than the processor elements ~ we expect the channel to 
be nearly independent of actuator decoherence. 
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FIG. 1. 2x(le-3n) Node Schematic. The nodes are taken 
to be identical, with resolved anisotropic hyperfine interac- 
tions (solid red lines) between electron actuator spins and 
nuclear processor spins. The local processors are initially 
disjoint, but may be effectively coupled (dotted lines) by 
modulating an isotropic actuator exchange interaction (solid 
blue double line) and moving into an appropriate microwave 
Hamiltonian interaction frame. The spin labeling is e^ for 
electron actuator spins and riij for nuclear processor spins, 
where i labels the nodes and j labels the qubits. 



The details of the method will be discussed from the 
standpoint of a spin-based, distributed node quantum in- 
formation processor. This system contains all the neces- 
sary physics and is representative of many other modal- 
ities being considered for experimental realizations of 
quantum information processing. For example, Rydberg 
atom excitations of neutral atoms 10 ;_ 11| , inductive cou- 



12, |13| , and Bloch wave 
all take the form of 



pling of superconducting qubits 

dispersion in cavity devices IJ, 

an isotropic dipolar coupling. Direct dipolar interactions 

also naturally occur in spin-based devices such as semi 

conductor quantum dots 



m 17 



^ .^, silicon-based devices 

la [iSJ , nitrogen- vacancy defect centers in diamond [19| , 
and other solid-state spin systems 2Cll422l|. The methods 
developed in this work may be readily extended to these 
systems. 
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FIG. 2. Energy Structure of 2x(le-ln) System. Quan- 
tum information is encoded in the states j4,4,00), |44-01), 
m 10), and m 11) of the actuator ground-state (Compu- 
tational) manifold, where arrows indicate electron actuator 
spin states and binary digits indicate nuclear processor spin 
states. The desired transition (bold red arrow) for a SWAP 
operation is implemented by applying a selective microwave 
field that induces transitions between the two manifolds (dot- 
ted lines) , efi'ectively moving the induced cross-node processor 
transitions (dashed lines) in the actuator zero-quantum (ZQ) 
manifold to the Computational manifold. Note that the ac- 
tuator excited state manifold, \tt), is not included as it is not 
involved in the internode transfer process. 



In our spin-based model, each node consists of a single 
actuator electron spin coupled via resolved anisotropic 
hyperfine interactions to each of k qubits of a local nu- 
clear spin processor (Fig. [1]) . Control over the local pro- 
cessors is achieved via electron-only modulation [23| , tak- 
ing advantage of the relative strength of the hyperfine in- 
teraction to generate a universal set of fast quantum gates 
on the nuclear spins 2J| , which serve as excellent storage 
elements for quantum information due to their relatively 
long coherence times [25| . The internode coupling of ac- 
tuators is given by an isotropic dipolar or exchange inter- 
action between electrons. The spatial separation of the 
nodes is taken to be sufficient for any cross-node dipolar 
interactions of nuclear spins to be negligible. 

The state structure of a two node system with one elec- 
tron actuator spin and one nuclear processor spin each 
- a 2x(le-ln) system - is shown in Fig. [2] The com- 
putational basis states of the local quantum processors 
are defined in the ground-state manifold of the actuators. 
This choice of encoding allows us to implement gates be- 
tween the disjoint processors by taking advantage of an 
induced cross-node coupling in the zero-quantum (ZQ) 
manifold of actuator excited states not used for informa- 
tion storage. 

We derive the form of the cross-node coupling for a gen- 
eral 2x(le-fcn) system. The nodes are taken to be iden- 
tical with energy structure given by a dominant, quan- 
tizing electron Zeeman interaction, H^, with a strong 
static magnetic field oriented along the laboratory z di- 
rection; a corresponding nuclear Zeeman interaction, H^; 



an anisotropic hyperfine interaction between electron and 
nuclear spins, H^^; and a dipolar interaction between 
electron spins, i?Q°- 

fj2c-kn = ^c ^ mi + H, 



rc-n j^ zjc-c 



(1) 



In a frame rotating at the electron Zeeman frequency, the 
resulting secular Hamiltonians are given in terms of the 
usual spin-i Pauli operators as [26| 



i?|=^w,^«-+a--) 

k 

i/S' = ^#.(a^^a"-+a^^a"-), 

k 



(2) 



where the vectors A^ = A^x + A^y + A'^z represent the 
strengths and directions of the hyperfine coupling be- 
tween the A:th nuclear spin in each node and the cor- 
responding actuator, w^' is the strength of the nuclear 
Zeeman interaction for the fcth nuclear spin, a;^ is the 
strength of the dipolar interaction, and a = (JxX + Uyij + 

<JzZ. 

The nuclear spins are quantized in an effective field 
given by the vector sum of the hyperfine and nuclear 
Zeeman interactions. The resulting eigenstates are non- 
commuting, allowing for universal control of the nuclear 
spins via electron-only control 2J|. Note that, since the 
nodes are identical, the above Hamiltonians do not pro- 
vide the ability to selectively address nodes. Universal 
control over the entire 2x(le-fcn) system is obtained by 
adding a term to the Hamiltonian that spatially labels 
the nodes to allow for local operations. These terms are 
not included in the present discussion as they are not 
necessary for the implementation of gates between nodes, 
and may be effectively turned off. We only require that 
the differences in the hyperfine coupling strengths within 
each node are large enough for each pair of identical spins 
to be spectroscopically resolved. This requirement lim- 
its the number of qubits per node [8| but, due to the 
inherent inefficiency of designing control sequences for a 
large number of particles, it is advantageous to keep the 
size of nodes small and rely on the ability to swap qubit 
states between nodes to implement large-scale quantum 
algorithms. 

The full set of interactions accessible by evolution un- 
der the Hamiltonians in ^ is given by the Lie algebra 
generated by taking Lie brackets to all orders 27J, |28j . In 
particular, the second-order bracket, [[-ff^ *', iJ^p], i/^p], 
takes the form of a four-body inter-node interaction, 
given by an effective cross-node nuclear spin dipolar cou- 
pling, ii/g~", along with fiip-flop transitions of the elec- 
tron spins: 



[[H^ ", H^p] , H^p] oc ojd 



( Ci Co 



^S""- (3) 



The resulting nuclear spin dynamics in the ZQ manifold 
may be decomposed into a sum of coupling terms, -ffp™, 





FIG. 3. Internode Coupling Network. The induced cou- 
pling network of local processor elements consists of interac- 
tions between every pair of cross-node spins. For implementa- 
tion of a parallel swap operation, we wish to keep interactions 
between identical spins (solid red lines) while refocusing all 
other interactions (dotted lines). 



that act on every pair of cross-node spins, nu and n2,n 
(Fig. [21). Each H^ may be written in terms of the 
well-known dipolar alphabet, with a± = a^ ^ iuy [26] : 
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(4d) 
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After application of an appropriate microwave control 
field, these four-body interactions in the ZQ manifold 
appear as effective cross-node two-body couplings in the 
computational manifold. To demonstrate this, we con- 
sider the implementation of a particularly powerful oper- 
ation: a parallel swap of entire local processor states be- 
tween nodes at once. For the case of a single nuclear spin 
per node, the relevant ZQ transitions are |ti 01) {]/[ 10| 
and Iti 10) (4,t 01|. Application of a microwave field with 
matrix elements |ti) (4.-H and \]/[) (4.4.I of strength com- 
mensurate with the ZQ transitions transforms both of 
the transitions to a swap operation in the computational 
manifold: lU 01) (U 10| (See Fig. [2). 

When multiple nuclear spins are present in each node, 
a parallel swap operation requires suppressing couplings 
between non-identical spins (^ 7^ m) while retaining cou- 
plings between identical spins (£ = m). This may be 
accomplished by exploiting the difference in symmetry 
between the prefactors of the coupling operators for iden- 
tical versus non- identical spins. 

Consider, as an example, a 2x(le-2n) system. The 
effective Hamiltonians of the induced interactions may be 



H^ + ij22 and H^ 
H^ . The effective dipolar coupling strength for 



written as H'^ + Hj^ , where TJq 

tt12 



H^ appears as odd order in A^ and A^, while each term 
in H^ appears as even order. Thus, by inverting the state 
of only the second (or first) spin in each node halfway 
through free evolution under the induced Hamiltonians, 
we can generate a zeroth-order average Hamiltonian of 
only the desired H'^ interactions [29| . Higher order terms 
in the average Hamiltonian may be suppressed through 
the use of more sophisticated pulses or by applying the 
evolution-pulse-evolution cycle at a rate fast compared to 
cjd- This symmetry argument may be easily generalized 
to a larger number of nuclear spins per node by applying 
a binomially expanding set of inversion pulses to properly 
select the desired couplings [30|. 

We now consider how to isolate the desired interac- 
tion ([3]) from other elements of the algebra. One method 
of suppressing the extraneous terms is to use a compos- 
ite pulse sequence to generate an effective Hamiltonian 
for which the desired second order commutator is the 
dominant term. Concretely, recall that by the BCH ex- 
pansion, e^e^ = exp(A: + Y+^[X,Y] + ^[X, [X, Y]] + 
rj^[F, [Y, X]] + •••). By recursively applying this expan- 
sion, we can derive an identity that suppresses all terms 
below second order: 

e^e^e-^e-^e-^e^e^e-^ = e[^'[^'^ll+-. (5) 

By making the correspondence X = ^hf ^^'^ ^ = -^d ° j 
we obtain a pulse composed of sequential periods of only 
electron dipolar or hyperfine evolution, leading to the 
effective propagator. 
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»r-'[[ffS"°,HS'].-HHF] 



(6) 



where higher-order terms have been neglected. We may 
also suppress the undesired terms by numerically op- 
timizing experimentally robust microwave pulses which 
achieve the desired interaction while suppressing all other 
interactions [3l|-|3J|. 

A final consideration is the sensitivity of the induced 
channel to actuator noise processes. We claim that by 
never transferring complete qubit state information to 
the actuators, we may operate in a regime where any 
portion of the information present in the ZQ manifold 
arrives back to the computational manifold before it is 
corrupted. We may quantitatively determine the robust- 
ness of the channel to actuator noise by comparing, as 
a function of noise strength, the Hilbert-Schmidt inner 

product fidelity, F(S'idcai, 5'noisy) = Tr {Sl^^^^S noisy) /d'^, 
between d? dimensional superoperators representing the 
channel in the presence of noise, S'noisy, and in the ideal 

noiseless case, 5'ideai- 

The ideal channel is generated by a Liouvillian oper- 
ator, L, corresponding to unitary evolution only. The 
noisy channel includes two dissipation operators, Di and 
D2, describing the relaxation of ei and 62, respectively: 



M 



-itL+tDx+tD2 



(7) 
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FIG. 4. Channel Noise Robustness. A plot of the channel 
fidelity, F, as a function of noise strength, Ti = T2, scaled by 
computing — log[l — _F]. The induced channel (solid line) per- 
forms significantly better than a serial swap operation (dotted 
line) [30]. Given a modest Rabi frequency, cji = 100 MHz, 
the induced channel is nearly independent of actuator deco- 
herence for electron relaxation times above 100 /is. 



A physically motivated model of noise is a contribution 
of phase and amplitude damping applied seperately to 
each electron, which leads to a dissipator, 



D-liT. 



Ta) (S_ 



Tia^ 






(8) 



where E^ = |0) (0| and E_ = |1) (1| are projection op- 
erators. The noise strength is parameterized by Fi 2, 
which are related to the commonly used energy, Ti, 
and coherence, T2, relaxation times by Fi — 1/Ti and 
T2 = i2Ti - T2) I (T1T2). A plot of the superoperator 
fidelity versus noise strength is shown in Fig. U) The 
noise has a minimal effect on the operation of the chan- 
nel for values of TiUJi > 10^. Assuming a modest Rabi 
frequency of oji = 100 MHz, actuator relaxation times of 
100 fj,s are required to avoid significant corruption of the 
information during transfer. Currently achievable relax- 
ation times for electron spins are well within this range 

By taking advantage of the additional degrees of free- 
dom present in an actuator based system, and manipulat- 
ing the naturally occuring actuator interactions between 
nodes, we were able to create an effective channel be- 
tween initially disjoint local processors that allows the 
parallel transfer of fc-qubit states between nodes, effec- 
tively independent of actuator decoherence. The chan- 
nel takes the form of four-body cross-node interactions 
in the zero-quantum manifold of actuator states which, 
after moving into an appropriate microwave interaction 
frame, appear as effective two-body couplings of cross- 
node qubits in the computational manifold. The result- 
ing complete cross-node coupling network may be used 
to generate a universal set of gate operations between 
nodes. We expect the techniques described in this work 
to be applicable to a wide-variety of quantum devices. 



with minimal need for modification. 
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In this supplemental material, we elaborate on the calculations used in the main text to justiiy the robustness of 
our proposed wide quantum channel against relaxation processes acting on the electron actuators. In particular, we 
compare the pulse sequence for isolating the second-order commutator of two operators given in Equation (5) of the 
main text to a sequence of pulses implementing serial swap gates. For simplicity, we shall assume that there is one 
nuclear spin in each node, and shall hence drop the node index from our nuclear spin labels. 

Concretely, by alternately suppressing and inverting the hyperfine and dipolar couplings, we can implement a 
unitary U generated by the effective cross-node nuclear dipolar coupling for a cycle of length tc — 8t, 

U{8t) ^ e^'^""^^ e^'''"'^'^ e+""'S^ e+""'°'^ e^'^"'S^ e^'''"'°'^ e^'''"'^^ e^""'^'^ 



e 

„3 



HF • -^-'HF '-"-^D 



= exp (+ir3wd«V!_^ + cr'L'cTf) ® i/g"") . 

A reasonable choice for tc, and hence for r, is such that the nuclear spins evolve for a phase tt/2 under the action 
of the effective nuclear dipolar interaction. We can estimate this by choosing r such that 

^3 II Tjc-e II II Tje-n ||^ ^ /o 

^ ||-"d II II-'^hfII - 1^/^- 

Since we are only interested in how this pulse sequence acts on the states of the nuclear spins, we compose the 
unitary evolution with superoperators representing state preparation Sprepip"^"^] = |ieiie2) {ieiie2\'^P^^^^ and partial 
tracing 5trace[p''^''^"^'^^] = Treie2 (p''i''2nin2^ q£ |.]-^g clcctrou spins. We also include a selective microwave Hamiltonian 
to cause mixing between the computational and zero quantum manifolds, as illustrated in Figure 2 of the main text: 

where E^^ E"_^ = |01) (01 1 and E"_^ E^ = |10) (10| are projection operators on the nuclear spins. The action of the 
ideal wide quantum channel is given by 

C ., , — Q . I p-*'''-£'HF+*''"-'^pw „-JrL53"+iTLpw p+irLgp+JrLpw p+irLp "+iTLjiw 

^woc. inna '-'trace I ^^ ^ '-■ ^ 



wqc, ideal *^ trace 

+iTLg"p+iTLj,w -iTL^"+iTLj,w ^-irL'^^+iTL-^^ „+iTL^"+iTLj,w ) c 
c; c; D D I ^prep; 

where Ly>, iuF and ipw are the Liouvillian representations of the respective Hamiltonians _ffD, -ffuF and -ffpw. In the 
column-stacking basis for Liouville space, we can write that 

L^t®H -H'^ ®t^t®H -H* ®t, 

where H* is the complex conjugate of H, so that L[p\ ~ Hp — pH = [H,p] represents the Liouville- von Neumann 
equation dtp ~ — iL[/C)]. 

To model how decoherence affects this sequence, we act on the electrons with the dissipator given in Equation (8) 
of the main body. This dissipator is then added to the generator for each interval in the above pulse sequence, so that 

C _ C [ —irL'^p+iTLy^+TD —irL'^'^+iTLy^^+Tb +irL'^p+iTL^,^^+Tb +iTLo "+irLy„+r_D 

^wqc, noisy '-^ trace I ^ ^ ^ ^ 

+iTL'^p+iTL^,^^+rD -irL'^"+iTL^,^+TD -iTL'^^+irL-^^+Tb +1x1,^ "+irLyw+T-D\ A 
o c o c I • Oprep. 

We are then concerned with how robust the action of our wide quantum channel is against the relaxation; that 
is, we are concerned with how accurately we may model our channel on the nuclear spins as decoherence-free. To 
answer this, we consider the process fidelity given by the Hilbert-Schmidt inner product between the ideal and noisy 

superoperators, F(5'wqc, ideal, S'wqe, noisy) = Tr(S'^q^ idear'5'wqe, noisy)/c?^, whcrc d is the dimension of the Hilbert space 
of nuclear spin states. 

For comparison to a serial swap protocol, we consider a composition of serial swap gates 

SWAP"!"^ = (SWAP^i'^iSWAP"^"^) • SWAP'=i'=^ • (SWAP"i°iSWAP"^°^) , 
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FIG. 5. Channel Noise Robustness. A plot of the channel fidelity _F as a function of noise strength, T\ = T2, scaled by 
computing — log[l — F]. The solid curve shows the fidelity of the wide quantum channel, while the dotted curve shows the 
fidelity of a serial SWAP protocol. Values of Hamiltonian parameters used for simulation are included in the table inset. For 
electron relaxation times above 100 ps, the channel is nearly independent of actuator decoherence. 



where we have grouped in parentheses those swap gates which can be performed in parallel. By noting that for two 
spins, SWAP = exp [f (^i?^ • f?^ — 1^1^)] , where ct' is the vector of Pauli matrices on spin i, we can represent SWAP"i"^ 
by the pulse sequence 



SWAp"i"2 =exp I thfOJM 
exp f tdUJd 
exp f thfUJhf 



.^ei . -ni _ tci-fim 



jcijr 



.-01 .^02 _icnc 



.r;C2 . ^.112 _ 1ie2tn2 



jeaji, 



where oJhf — H-f^HF lli ^^'^ where ihf and td are chosen to ensure that thfWhf — tdUJd — 7r/2. As before, we can represent 
this sequence as an action on the nuclear spins alone by representing each Hamiltonian as a Liouvillian, calculating the 
superoperators and composing with the preparation and trace steps S'prcp and 5traco- We model decoherence in the 
same way as before, adding the dissipator D to each step of this serial protocol. Doing so, we obtain superoperators 
•SswAP, ideal and S'swAP, noisy for the ideal and noisy serial SWAP protocols, respectively. 

In Figure 4 of the main text, reproduced here as Figure [SJ we compare the fidelities F(S'wqc, ideal, S'wqc, noisy) and 



F{Ss 



Dsv 



SWAP, ideal, »JSWAP, noisy 

decoherence than the serial swap protocol. 



^), observing that the wide quantum channel performs significantly better against electron 



